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ABSTRACT 

We present new results from our search for Giant Hii Regions in galaxies visible 
from the southern hemisphere. In this work we study two galaxies: NGC 7479 and 
NGC 6070. Using high-resolution spectra, obtained with different instruments at Las 
Campanas Observatory, we are able to resolve the emission-line profile widths and de- 
termine the intrinsic velocity dispersion of the ionised gas. We detect profile widths cor- 
responding to supersonic velocity dispersions in the six observed H ii regions. We find 
that all of them show at least two distinct kinematical components: a relatively narrow 
feature (between 11 and 22 kms~^) and a broader (between 31 and 77 kms~^) 
component. Two of the regions show a complex narrow profile in all ion lines, which 
can be further split into two components with different radial velocities. Whereas the 
wing broadening of the overall profile can be fitted with a low-intensity broad com- 
ponent for almost all profiles, in one region it was better reproduced by two separate 
shell-like wings. We have analysed the impact that the presence of multiple compo- 
nents has on the location of the Hll regions in the \og(L) — log((7) plane. Although the 
overall distribution confirms the presence of a regression, the precise location of the 
regions in the plane is strongly dependent on the components derived from the profile 
fitting. 

Key words: (ISM:) Hii Regions - galaxies: starburst - galaxies: individual: NGC 
7479 - galaxies: individual: NGC 6070 



1 INTRODUCTION 

The most extended and luminous objects which are observed 
on the disks of spirals, in irregulars and in starburst galax- 
ies are called Giant Extragalactic Hll Regions (GHliRs). 
GHliRs are places of very active star formation and they 
are characterised by the large emission of ultraviolet photons 
from a large number of young and massive stars. Owing to 
the large number of ionising photons most of the surround- 
ing gas is ionised, and the recombination process provides a 
strong signature of recent or ongoing massive star formation. 

The emission-line profile widths of the giant Hll re- 
gions imply the existence of supersonic motions in the gas 
(Smith & Weedman 1970). Furthermore, Terlevich & Meb 
nick (1981) found a correlation between the gas velocity dis- 
persion ((j) and the total luminosity (L) emitted in the re- 
spective line. Hippelein (1986) and Arsenault & Roy (1988) 
confirmed the existence of such a regression, but there was 
no agreement on the values derived for its slope so no conclu- 
sion could be reached regarding the origins of this motion. 
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Fuentes-Masip et al. (2000) and Bosch et al. (2002) using 
different techniques on different samples agreed on a slope 
close to a value of 4, which seemed to favour the gravitational 
source of energy for the observed supersonic motion. Even if 
the nature of this behaviour is not fully understood, it still 
allows us to distinguish giant Hll regions from an agglom- 
eration of classical Hll regions by means of high resolution 
spectroscopy. GHllRs provide the link between small-scale 
star forming regions, such as Orion in our Galaxy and violent 
star formation processes taking place in starburst galaxies. 

In Firpo et al. (2005) we have confirmed the giant na- 
ture of three candidates to giant Hll regions in the southern 
sky which were identified as very luminous Hll regions by 
Feinstein (1997). 

Here we continue our search and detailed analysis of 
GHllRs in local Universe galaxies. We have selected from 
the Feinstein catalogue, a sample of the brightest Hll regions 
from two spiral galaxies: NGC 7479 and NGC 6070. 

NGC 7479 is a barred spiral galaxy, at a distance of 
31.92 Mpc (Rozas et al. 1999), classified as SB(s)c by de 
Vaucouleurs et al. (1991). Devereux (1989) and Martinet 
& Friedli (1997) have classified it also as a starburst galaxy, 
whereas Keel (1983) and Ho et al. (1997) have identified it as 
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a LINER and a Seyfert 1.9 galaxy respectively. Many studies 
have been performed on this galaxy owing to the presence 
of a certain asymmetry in the spiral structure (Quillen et al. 
1995) that leads several authors to suggest this might be 
owed to a recent merger event (Laine 1996; Laine & Heller 
1999). This scenario can also account for perturbations in 
the velocity field (Laine & Gottesman 1998) and the lat- 
est minor merger model indicates that the remnant may 
be situated within the bar (Laine 2001) as no remnant of 
the merger could be identified in optical images (Saraiva & 
Benedict 2003). 

NGC6070, at a distance of 29.8 Mpc (Garcia-Gomez 
et al. 2002), is classified as SA(s)cd by de Vaucouleurs 
et al. (1991) although Grosb0l et al. (2004) state that a 
bar with two main arms can be identified. Pohlen & Tru- 
jillo (2006) further studied the bar structure classifying it as 
Type II.o.CT according to Erwin et al. (2008). The galaxy 
shows a normal arm distribution, with two arms that break 
into several on each side of the central region (see Figure 2). 
The disk dynamics was studied by Marquez et al. (2002), 
and Garcia-Gomez et al. (2002) studied the distribution of 
Hll regions in the galaxy. 

In order to verify the nature of our giant Hll region can- 
didates, we obtained high resolution spectra to measure the 
emission-line profile widths and to estimate if the velocity 
dispersion is indeed supersonic. To determine the optimum 
number of Gaussian components to fit each line profile in or- 
der to minimise the uncertainty in the velocity dispersion of 
the ionised gas, we have considered information from other 
studies available in the literature. Previous papers have dealt 
with the fact that a single Gaussian profile might not be 
realistic enough when fitting the observed emission lines. 
Several studies favour the presence of an omnipresent broad 
component which explains the integral profile wings. They 
include Melnick et al. (1999) who identified a broad compo- 
nent in the recombination lines in the central region of the 
30 Doradus nebula, Diaz et al. (1987) and Terlevich et al. 
(1996) in the M33 giant Hll region NGC604, Mendez & Es- 
teban (1997) in four Wolf-Rayet galajcies, and Homeier & 
Gallagher (1999) in the starburst galaxy NGC 7673. More 
recently, Hagele et al. (2007); Hagele (2008); Hagele et al. 
(2009) and Hagele et al. (2010) found that the best Gaus- 
sian fits in circumnuclear star-forming regions involved the 
presence of broad and narrow components for the emission 
lines of the ionised gas. 

A variety of line broadening mechanisms have been pro- 
posed to interpret the existence of the broad supersonic 
component measured in the emission line profile of GHllRs. 
Dynamics of virialised systems (Tenorio-Tagle et al. 1993), 
superposition of multiple gas bubbles in expansion (Chu & 
Kennicutt 1994), or turbulence of the same interstellar gas 
(Medina Tanco et al. 1997). More recently, Westmoquette 
et al. 2007; 2007 concluded that the narrow component rep- 
resents the general disturbed ionised interstellar medium 
(ISM), arising through a convolution of the stirring effects of 
the starburst and gravitational virial motions. On the other 
hand, the broad component results from the highly turbu- 
lent velocity field associated with the interaction of the hot 
phase of the ISM with cooler gas knots, setting up turbulent 
mixing layers. 

Other studies support the existence of shell-like profiles 
on the blue and red wings of the main feature as found by 



Chu & Kennicutt (1994), who state that the integrated line 
profiles of the entire 30 Dor region may be fitted with a 
broad Gaussian with low-intensity wings. Relano & Beck- 
man (2005) claim that a large fraction of the Hll regions 
in three barred spiral galaxies shows high- velocity and low- 
intensity features in its integrated line profiles and Rozas 
et al. (2006) find that the great majority of their sampled 
Hll regions show that the brightest principal component 
is accompanied by low-intensity components symmetrically 
shifted in velocity. 

In this paper we present echelle data obtained at Las 
Campanas Observatory (LCO) of six GHllRs candidates in 
the two mentioned galaxies. We have determined the ve- 
locity dispersion and we have evaluated the possible pres- 
ence of a broad component or two symmetric low-intensity 
components in the fit with the observed emission line pro- 
file widths. In section 2 we present the observations and 
the data reduction. Section 3 details the analysis performed 
on the emission-line profiles and discusses each conspicuous 
case individually. The summary and conclusions of this work 
are included in Section 4. 



2 OBSERVATIONS AND REDUCTIONS 

According to canonical denomination for GHllRs we have la- 
belled our candidates following their brightness ranking from 
Feinstein catalogue. NGC 7479 I is therefore the brightest re- 
gion in NGC 7479, NGC 7479 II is the second brightest and 
so on. Observed regions from NGC 7479 and NGC 6070 are 
shown in Figures 1 and 2 respectively. High resolution spec- 
tra were obtained with an echelle spectrograph at the 100- 
inch du Pont Telescope, Las Campanas Observatory (LCO) 
between 19 and 22, July 2006. The spectral range covered 
by the observations was from 3800 to 9500A. Observing con- 
ditions were good with 1 arcsec seeing and photometric sky. 
A 2x2 binning was applied to the CCD in order to min- 
imise readout contribution to the final spectrum noise. The 
spectral resolution achieved in our du Pont Echelle data 
(obtained with 1 arcsec effective slit width) is R~25000: 
AA=0.25A at A 6OOOA, as measured from the FWHM of 
the ThAr comparison lines. This translates in a resolution 
of .-12 kms"\ 

The standard stars, HR7950, and HR4468 and the 
CALSPEC spectrophotometric standard star, FeigellO 
(Bohlin et al. 2001), were also observed for fiux calibration 
purposes. In addition, Th-Ar comparison spectra, milky fiats 
(sky fiats obtained with a diffuser , during the afternoon) 
and bias frames were taken every night. The exposure time 
for standard stars was five seconds for the bright HR stars 
and 1200 seconds for FeigellO. 

We also obtained high resolution spectra for the three 
Hll regions in the NGC 6070 galaxy using the echelle double 
spectrograph Magellan Inamori Kyocera Echelle (MIKE) at 
the 6.5-m, Magellan II (Clay) Telescope, LCO, in July 2004. 
No binning was applied to the 2Kx4K CCD detector and a 
1 arcsec slit was used. The spectral resolution measured on 
our MIKE spectra is 11 kms~^, which is quite similar than 
that obtained at du Pont. The exposure time for CALSPEC 
spectrophotometric standard star BD-I-28D4211 was 900 
seconds. Th-Ar comparison lamps, milky flats (internal fiat 
field lamp with a diffuser glass slide) and bias frames were 
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Figure 1. Ha image witli the GHllR candidates observed in 
NGC 7479, identified by circles. Ha images were kindly made 
available by C. Feinstein. 
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used to calibrate these data. Table 1 lists the journal of 
observations, together with the exposure times, air masses 
and the standard stars that were used for flux calibrating 
each region. 

The data analysis was carried out with IRAF^ software. 
After bias subtraction and fiat field corrections with Milky 
Flats, the bidimensional images were corrected for cosmic 



^ Image Reduction and Analysis Facility, distributed by NOAO, 
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Table 1. Journal of observations for the candidate GHnRs and 
flux standard stars. The candidate regions have labelled following 
their brightness ranking from Feinstein catalogue (see Feinstein 
(1997)). 



Telescope 


Region 


Date 


Exp. 


sec z 


Standard 


du Pont 


74791 


2006 Jul 19 


1800 


1.6 


Feige 110 




7479 II 


2006 Jul 20 


1800 


1.5 


HR7950 




7479 III 


2006 Jul 20 


1800 


1.4 


HR7950 




60701 


2006 Jul 21 


1800 


1.2 


HR4468 




6070 II 


2006 Jul 21 


1800 


1.2 


HR4468 




6070 IV 


2006 Jul 21 


1800 


1.2 


HR4468 


Clay 


60701 


2004 Jul 12 


1800 


1.3 


BD+28D4211 




6070 II 


2004 Jul 12 


1800 


1.2 


BD+28D4211 




6070 IV 


2004 Jul 12 


1200 


1.7 


BD+28D4211 



rays with the task cosmicrays which detects and removes 
cosmic rays using a flux ratio algorithm. The corrected data 
were reduced by IRAF routines following similar procedures 
to those described in Firpo et al. (2005). 

We also compared the red end of the wavelength cali- 
brated spectra with the night-sky spectrum by Osterbrock 
et al. (1996). This turned out to be a very reliable confir- 
mation of the goodness of the wavelength solution, and we 
were able to check that differences between our wavelengths 
and the sky line wavelengths were below 0.05A. 

We performed flux calibrations using the observed spec- 
trophotometric standard stars as described in Firpo et al. 
(2005). The observed standard stars had their fluxes tabu- 
lated every 16A, and the amount of defined intervals within 
an echelle order ranged from four to twelve, depending on 
the quality of the spectrum. 

For the Magellan spectra, and in the case of the 
NGC 7479 1 region du Pont data, we were able to obtain 
CALSPEC spectrophotometric standard stars (Bohlin et 
al. 2001), whose fluxes were tabulated every 2A. These stars 
are ideal for calibrating high resolution echelle spectra, al- 
though their relatively low brightness (Feige 110, V=11.83 
and BD-f28D4211, V=10.51) makes them time-consuming 
targets for two metre class telescopes. 



3 ANALYSIS OF LINE PROFILES 

Making use of the known redshifts z=0.Q06685 for NGC 6070 
(Springob et al. 2005) and z=0.007942 for NGC 7479 
(Haynes et al. 1998) we identify the hydrogen recombina- 
tion lines, such as Ha and H/3, and coUisionally excited lines, 
such as [Nll]A6548,6584A, [Sll]AA6717,673lA, present in the 
spectra. These strong lines are used to analyse the structure 
of the profiles as they allow us to verify the existence of more 
than one component. The adopted laboratory wavelengths 
are taken from work by Garcia- Rojas et al. (2005). 

To determine the radial velocities and the velocity dis- 
persions of the ionised gas, we measure the central wave- 
length and width of several emission lines. The radial veloc- 
ity of each emission line is obtained from its central wave- 
length determined from the Gaussian profile fitting, and 
their errors are given by the fitting error provided by the re- 
sampling done within the NGAUSS task of IR AF and taking 
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into account the rms of the wavelength cahbration. The true 
velocity dispersion (a), of each emission line is calculated as 
described in Firpo et al. (2005). The width of the observed 
profile (o"o) is affected by the contribution of thermal random 
motions (at) and the instrumental profile (ci). Then, we can 
obtain the true velocity dispersion as a^—a^-af-at- As in 
Firpo et al. (2005), a typical kinetic temperature T = lO^K 
is assumed, and the instrumental profile (ai) is very well 
approximated by a single Gaussian function. The velocity 
dispersion errors have been calculated using the observa- 
tional errors in the observed profile and assuming negligible 
errors in ai and at- Owing to the generally high metallicity of 
these kinds of objects as a class, the temperatures are prob- 
ably lower than lO^K (see Di'az et al. 2007, and references 
therein). The low excitation of these objects estimated from 
their spectra makes any temperature sensitive line too weak 
to be observed. In fact, in most cases, the [Oin]A5007A fine, 
which is typically one hundred times more intense than the 
auroral [Oni]A4363A one, can barely be seen. In this way, 
we have also calculated the velocity dispersions assuming a 
kinetic temperature of 5000K and we do not find great dif- 
ferences with the previous calculations derived using 10*K. 
These differences are between 0.02 and 3 kms~^ depending 
on the measured velocity dispersion since (Jt is constant for 
a given temperature and subtracted in quadrature to CTq. 
Obviously, the estimated a are always higher for the lower 
assumed kinetic temperatures. Then, for consistency with 
Firpo et al. (2005) we use a kinetic temperature of 10'*K. 
Therefore, this effect owing to the kinetic temperature un- 
certainty, would produce an underestimation in the derived 
true velocity dispersions. 

When fitting single Gaussian profiles to the emission 
lines we notice that, although it works well for most of the 
observed profiles, there is a residual present in the wings 
of several lines. This has been observed in GHllRs before. 
Some authors (Melnick et al. 1999; Hagele et al. 2007, 2009, 
2010, among others) have shown that the profile wings are 
well fitted by single (broad) Gaussian of supersonic widths, 
whereas other authors (Chu & Kennicutt 1994; Relafio & 
Beckman 2005; Rozas et al. 2006) proposed that the profiles 
appear to be dominated by a few, well-resolved, expanding 
shells. Wherever possible, we have considered and applied 
both suggestions. 

In order to fit multiple components to the observed 
profiles we make use of the ngaussfit task within IRAF. 
Ngaussfit performs iterative fitting of multiple Gaussian 
profiles to spectral lines. The task allows us to select indi- 
vidual parameters to be fitted and delivers an estimate for 
the uncertainty of derived parameters. The task needs initial 
guesses for the function coefficients, and these are specified 
through previously generated tables. 

In those cases when more than one component is evi- 
dent, we follow an iterative procedure in which we only allow 
the task to fit a limited subset of parameters at one time. 
This is done to constrain the universe of possible solutions 
to the fit by making use of the different information avail- 
able from each emission line. Basically, after obtaining an 
initial guess for the second component parameters, we fit 
the amplitude, centre and width of each component at one 
time, leaving the parameters of the other one untouched. 
Once the central wavelength of all components is relatively 
well known, we set them as fixed parameters and allow the 



task to perform a final fit of amplitude and width of all 
components at the same time. 

The profile fittings show the presence of a residual in 
the emission line wings for all cases, becoming more evi- 
dent in the strong Ha emission line. Therefore, we evaluate 
the presence of either a broad component or shell features. 
For the former assumption, we introduce a broad compo- 
nent with an initial width which is three times the width of 
the narrow component (s). This guess comes from an over- 
all view of the values found by Hagele et al. (2007) when 
performing a similar analysis. We fit the peak values for the 
narrow components and then their profile widths. The broad 
component parameters are then fitted before we start the it- 
erative procedure again. In the latter scenario, a similar it- 
erative procedure takes place, but in this case we add a blue 
and a red narrow component to account for the additional 
emission. The validity of the profile multiplicity and broad- 
ening is checked over the different emission lines available 
for each region, although the most reliable profile analysis 
comes from the strongest lines, such as Hq. 

Out of our sample of six GHllRs we find that all of 
them show evidence of wing broadening. For these regions 
we attempt to fit an overall broad component or two narrow 
components symmetrically shifted in velocity with respect to 
the intense component. In the following subsections we will 
discuss our findings grouping them by the outcome profiles. 

3.1 NGC7479I 

Rozas et al. (1999) undertook a complete study to date of 
the Hll regions in NGC 7479 galaxy, determining positions, 
angular sizes, and absolute fiuxes of over 1000 Hll regions 
and they constructed the luminosity function for the regions 
over the whole galaxy. We have cross-correlated these results 
with our sample and find that the three NGC 7479 regions 
selected in our work can be identified in the Rozas cata- 
logue: NGC7479I as RZH99=1, NGC7479II as RZH99=7 
and NGC7479 HI as RZH99=74 (RZH99 is the number iden- 
tification used by Rozas et al. 1999). 

We identify and fit the Gaussian profiles to the Ha, 
[Nn]A6584A and [Sii]AA6717,673lA lines in NGC 74791 re- 
gion with splot routines. The Gaussian fits reveal the pres- 
ence of two distinctly separated kinematic components (la- 
belled A and B) in all analysed profiles. Component A shows 
a much broader profile than component B, and both compo- 
nents show little spread in their individual radial velocities 
among different emission lines present in the spectrum. The 
reliability of these values is confirmed when we improve the 
profile fitting using NGAUSSFIT which in turn yields values 
for the profile width of each component. 

Even considering the presence of two distinct compo- 
nents with different radial velocities, the overall fit cannot 
account for a residual emission in the profile wings. Fol- 
lowing the procedures outlined in the previous section, we 
are able to fit a broad component, with a velocity disper- 
sion of about 75 kms^^ from the Ha emission line and 
about 60 kms^^ from the other ones. An attempt to fit 
twin narrow profiles did not converge to a reliable solution. 
Table 2 shows parameters for the three components that fit 
the global profile. Individual Gaussian component fluxes are 
listed as emission measure (EM) relative to the total line 
flux following the work by Relafio & Beckman (2005). The 
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Figure 4. 2D Ho velocity image of NGC7479I region with the 
[Nll]A6584A (red) contours overlapped. Three sections are num- 
bered in the 2D image. 



overall Ha flux, uncorrected for reddening, is found to be 
3.37 X 10"^"' ergs^^cm"^. In Figure 3 we show the NGAUSS 
fitting done with three different Gaussian components in the 
emission lines which have enough signal to provide an accu- 
rate fit. The emission lines with a low signal-to-noise ratio 
for which the fits do not provide accurate results for the 
broad component, are not listed in the table. The velocity 
dispersion for the ionised gas is derived by taking Oi — 5.1 
kms~^ as instrumental width, and a kinetic temperature of 
T ~ W^K as mentioned above. 

Since part of the spatial information along the slit is 
lost through the aperture extraction process we generate 2D 
velocity images in Ha, [Nll]A6584A and [Sll]A6717A lines. 
Taking into account that the slit size is l"x 4", the data 
are binned by two pixels in both directions and the seeing 
was 1", we take aperture sizes of 0.5 arcsec ~ Ipx, there- 
fore dividing the slit into eight sections. We extract traces 
of the spectrum at eight representative cuts along the slit, 
preserving the same spatial resolution. This is done for each 
echelle order where a line of interest is present and each 
spectral section is rebinned to velocity space. In this ref- 
erence frame, individual slices are stacked together using 
the imstack task, which allows us to generate 2D velocity 
images. Figure 4 shows the 2D Ha velocity image with the 
[Nn]A6584A contours overlapped. It can be readily seen that 
individual knots can only be identified when the kinemati- 
cal information is considered in this Vr vs. slit section plane. 
This suggests that 3D spectroscopy will be able to provide a 
powerful insight on the structure of this type of star forming 
knots. 

We have performed a detailed analysis on the slit stack, 
by fitting Gaussian profiles at each slit section. The result- 
ing parameters are listed in Table 3 and the graphics output 
can be seen in Figure 5 where is shown 2, 4, 5, and 7 Ha 
sections, and 4 and 5 [Nll]A6584A sections. Focusing on the 
slit edges first (sections 1 to 3 and 6 to 8) we identify two 
distinct kinematical components at 2500 kms~^ and 2540 
kms~^ labelled A and B respectively. Each of these rela- 
tively narrow features has a broader component associated 
with it. Although the narrow component can still be de- 
blended in the central region, we are not able to disentangle 
the broad ones, which results in a single broad profile with 
a velocity similar to the average between components A and 
B. Regarding the velocity dispersion of each component, we 
have listed them in Table 3 as a reference, but they should 
not be taken as a definite value as the signal-to-noise ratio 
is variable and somewhat low for the slit edges. 

If we compare these results with the radial velocity map 
derived for the ionised gas in NGC 7479 (Laine 1996) we 



Table 3. Gaussian fitting results for individual sections of the 
Hq emission line profile in NGC 7479 I. Kinematical parameters 
(radial velocity (Vr) with errors ( kms~^), the intrinsic veloc- 
ity dispersions {(Jint) corrected for the instrumental and thermal 
widths ( kms~^), and emission measure (EM) of each component 
as a fraction of the total EM of the sampled section (%) data are 
shown. 



Ao 


section 


comp 




error 




error 


EM 


6563 


Ha-1 


A 


2501.8 


0.9 


12.6 


1.2 


51 






broad 


2500.7 


4.8 


59.3 


4.3 


49 


6563 


H«-2 


A 


2503.0 


0.9 


11.6 


1.5 


63 






broad 


2500.3 


1.6 


37.7 


2.5 


37 


6563 


H«-3 


A 


2504.7 


0.6 


16.4 


0.8 


43 






broad 


2506.8 


1.4 


47.0 


1.7 


57 


6563 


Ha-4 


A 


2500.8 


0.7 


15.8 


0.9 


43 






broad 


2519.6 


2.9 


55.0 


3.6 


30 






B 


2536.8 


0.7 


10.7 


0.9 


27 


6563 


Ha-5 


A 


2500.9 


0.6 


10.5 


0.8 


25 






broad 


2515.6 


1.5 


55.9 


1.9 


35 






B 


2535.6 


0.6 


13.1 


0.8 


40 


6563 


Ha-6 


broad 


2502.6 


1.5 


45.5 


1.7 


39 






B 


2538.6 


0.6 


17.4 


0.8 


61 


6563 


Ha- 7 


broad 


2519.4 


1.2 


46.0 


1.3 


56 






B 


2540.7 


0.6 


6.0 


0.8 


44 


6563 


Ha-8 


broad 


2523.9 


2.1 


44.9 


1.9 


80 






B 


2532.6 


0.6 




1.2 


20 



find, from visual inspection, that the expected velocity at 
the position of NGC7479 1 ranges between 2520 and 2535 
km s~ ^ close to an average value of our components A and B. 
A similar comparison, using the H I map by Laine & Gottes- 
man (1998) which has a better spatial and spectral resolu- 
tion, yields a value between 2530 and 2540 kms~^ much 
closer to the value measured for component B. Although the 
latter map is sampling neutral rather than ionised hydrogen, 
the former might be affected by the complex structure we 
find within the region. At this stage we can only suggest 
that component A shows an odd kinematic behaviour and 
that a detailed study of NGC7479 1 and its surroundings is 
needed to clarify this issue. 



3.2 NGC7479II 

A similar procedure is used to analyse the profiles of the 
emission lines of NGC7479 II. In this case, two compo- 
nents are identified with different radial velocities, but as 
the overall signal obtained for this region is lower than for 
NGC7479 1 it is not possible to disentangle them with 2D 
spectra. These two narrow components are found in the 
Gaussian profiles of the Ha, H/3, [Nil] AA6548,6584A and 
[Sn]AA6717,673lA lines. In this region, however, the wings of 
the brightest component profile could not be fitted with one 
single broad component. On the other hand, we are able to 
fit it with two wings only for the two lines with the strongest 
signal, namely Ha and [Nn]A6584A (see Figure 6). The fi- 
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Table 2. Results of Gaussian profiles fitting to the observed emission lines in NGC 7479. Each emission line is identified by its ion 
laboratory wavelength and ion name in columns 1 and 2. According to the different fits performed on each line, column 3 identifies each 
"narrow component" (A and B, where applicable), a broad component and/or blue and red wings (b wing and r wing respectively). Radial 
velocities (Vr) and intrinsic velocity dispersions {(Tint) together with their respective errors arc expressed in kms~-^. The intrinsic velocity 
dispersions are corrected for the instrumental and thermal widths. Emission measures (EM) are included and shown as a percentage of 
the component flux relative to the total EM of the region. 



NGC 7479 I NGC 7479 11 NGC7479 IU 





ion 


conip. 


Vr 


error 


0"mt 


error 


EM 


Vr 


error 


Tint 


error 


EM 


Vr 


error 


Tint 


error 


EM 


6563 


Ha 


A 


2500.4 


0.7 


22.3 


0.9 


39 


2492.4 


0.6 


14.8 


0.8 


66 


2481.1 


0.7 


11.3 


2.9 


67 






r> 


2538.9 


0.6 


15.3 


0.8 


42 


2532.2 


0.7 


11.1 


0.9 


23 
















broad 




o.y 


/ D.D 
















9479 ^ 


9 1 
z. 1 


J .0 


n Q 
u.y 


00 






b wing 












2442.7 


1.7 


13.5 


2.0 


5 
















r wing 












2573.8 


1.9 


18.0 


1.7 


6 












6584 


[Nil] 


A 


2506.8 


0.7 


21.6 


0.9 


39 


2494.2 


1.1 


16.4 


1.0 


53 


2482.9 


0.8 


12.8 


2.3 


52 






B 


2540.1 


0.7 


15.2 


0.8 


39 


2532.5 


2.0 


16.7 


1.3 


30 
















broad 


2523.4 


4.4 


63.6 


3.6 


22 












2478.4 


4.0 


40.4 


1.1 


48 






b wing 












2459.8 


9.1 


24.8 


4.3 


13 
















r wing 












2581.5 


6.5 


14.7 


2.5 


4 












6717 


[SII] 


A 
B 






















2480.9 


0.8 


10.2 


5.3 


33 






broad 






















2471.1 


4.1 


52.2 


1.4 


67 






b wing 




































r wing 
































6731 


[SIX] 


A 


2502.4 


1.1 


19.3 


1.2 


30 












2481.0 


0.8 


15.4 


4.3 


50 






B 


2538.0 


0.8 


14.3 


0.9 


39 


























broad 


2516.2 


4.2 


52.1 


5.1 


31 












2454.7 


6.4 


51.8 


1.2 


50 






b wing 




































r wing 

































Figure 3. Ngauss fits with three Gaussian components in the NGC7479I emission line profiles, two narrow principal components and 
a broad one. In order: Ha, [Nll]6584A, and [Sll]673lA. 
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nal fitting parameters are shown in Table 2. The overall 
Hq flux, uncorrected for reddening, is found to be 1.06 x 
10"^* ergs~^cm~^. The emission lines with a low signal-to- 
noise ratio for which the fits do not provide accurate results 
for the broad and the two wings components, are not listed 
in the table. 

We compare the velocity separation between the red 
(and blue) shifted components with the average of the cen- 
tral features and the Ha luminosity of the region. The veloc- 
ity separation of the red and blue shifted components in this 
work is in agreement with the regression found by Relafio & 
Beckman (2005) for their Hll regions sample. 

We also analyse the radial velocity behaviour within the 
galaxy, and observe a similar situation to that of NGC7479 1. 
There is no difference, within the errors, when we compare 
our average value derived from the Hq emission line with 
the Hq velocities field map (Laine 1996). If, however, the 
comparison is done against the HI rotation velocity map 



(Laine & Gottesman 1998), knot B seems to follow galactic 
rotation whereas knot A does not. 

3.3 NGC7479III 

The observed spectrum of this region shows, although weak, 
the [Olll]AA4959,5007A lines. This is probably owed to a dif- 
ference in metallicity, NGC7479 HI being slightly less metal- 
lic than the brightest two regions discussed above, or a dif- 
ference in the ionisation structure. The latter might be due 
to a relatively small difference in the hardness of the ionisa- 
tion radiation field (and, what is equivalent, in the ionising 
stellar population) responsible for the excitation observed 
in the spectra or a difference in the nebular geometry (see 
a detailed discussion in Diaz et al. 2007; Hagele et al. 2006, 
2008; Perez-Montero et al. 2010, and references therein). 
Unfortunately, as we said above, since no temperature sen- 
sitive lines can be observed, we are not able to derive their 
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metallicities, and distinguish the origin of this small diflter- 
ence in the excitation present in the spectra of these regions 
of NGC 7479. 

We identify and fit the Gaussian profiles to the 
[Oiii]A5007A, Ha, [Nn]A6584A and [Sn]AA6717,673lA lines 
in NGC7479III region with ngaussfit. In this region we 
have found one supersonic component, with no evidence of 
multiple contributions but it is possible to detect flux excess 
on the wings of the emission lines. We therefore fit a broad 
Gaussian component together with the narrow component 
(see Figure 7). This procedure works well for almost all de- 
tected lines. In the case of [Oiii]A5007A, however, there are 
no prominent residuals when we fit only one Gaussian com- 
ponent which has a dispersion of about 20 kms~^and we are 
therefore not able to fit a broad Gaussian component. For 
the other listed lines the derived radial velocities, the veloc- 
ity dispersions (using af—^A kms~^ and T ~ lO^X) and 
errors are listed in Table 2. The overall Ha fiux, uncorrected 
for reddening, is found to be 5.43 x 10~^^ ergs~^cm~'^. It 
is worth mentioning that there seems to be some systematic 
offset between the derived profile centres of the narrow and 
broad profiles which is, in average, about 10 kms^^. This 
happens for all lines, although it becomes more evident for 
[Sn]A673lA. An undiscovered blueshifted narrow component 
could account for the observed behaviour, although better 
signal-to-noise data are needed to confirm this hypothesis. 

The radial velocity average derived for the single com- 
ponent is in agreement with the expected radial velocity at 
the position of NGC7479 HI, showing only a slight offset of 
(10 kms~^), well within the uncertainties. 

3.4 NGC6070 I, II & IV 

All regions observed in NGC 6070 show one narrow com- 
ponent together with the underlying broad component. 
We identify and fit the Gaussian profiles to the H/3, 
[Oiii]A5007A, Ha, [Nn]A6584A and [Sii]AA6717,673lA lines 
in almost all regions, (see Figures 8-10). Table 4 shows the 
derived kinematical parameters (using erf — 5.3, 5.6 and 
5.3 kms~^ for regions I, II and IV respectively; T is set to 
~ lO^Js' for all regions). 

For these three regions, we are also able to compare 
the du Pont data with higher resolution spectra data ob- 
tained with the double echelle spectrograph MIKE in July 



2004. The agreement in the kinematical parameters can be 
seen in Table 4, which shows the reliability of the data here 
presented. The overall Ha flux, uncorrected for reddening, 
for du Pont data are 1.68 x 10~", 6.26 x 10"" and 2.15 
X lO"^'* ergs^^cm"^ for regions I, II and IV respectively. 
And for MIKE data are 1.74 x 10~^^ 1.29 x 10"" and 1.88 
X 10"^'^ ergs"^cm"^. Only a slight bias towards lower ve- 
locity dispersion can be seen for NGC6070I, which can be 
probably attributed to differences in the centering of the 
slit during acquisition of the target. There seems to be some 
systematic offset between the derived profile centres of the 
narrow and broad profiles. Theses values are, in average, 
about -25 kms"^ for du Pont data and -10 kms"^ for MIKE 
data in the case of NGC 6070 II, and about 10 kms"^ and 
5 kms"i for NGC 6070 IV. 

Regarding radial velocities of the three observed regions 
in NGC 6070, there is no detailed velocity map in the liter- 
ature. We are able to compare them with rotation curves 
derived from long slit observations made by Marquez et al. 
(2002). Although the orientation of the slits does not match 
exactly the location of all our regions, the velocities are in 
fair agreement considering the angular distance of the Hll 
regions to the centre of NGC 6070. 

3.5 Relation between Ha Luminosities and 
Velocity Dispersion 

We have also analysed the impact that the presence of mul- 
tiple components has on the location of Hll regions in the 
log(I/) — log(cr) plane. To evaluate this we have plotted 
in Figure 11 the luminosities and velocity dispersions de- 
rived directly from our spectrophotometric data. Luminosi- 
ties were derived from the fluxes measured directly from 
the component fltting to our echelle spectra (uncorrected by 
reddening), and using distances as published by Rozas et al. 
(1999) for NGC 7479 and by Garci'a-Gomez et al. (2002) for 
NGC 6070. We are aware that slit spectra might not provide 
the most realistic sample of the overall flux of the ionised 
nebula, but the accuracy of acquisition and guiding is not 
good enough as to allow us to perform a proper identifi- 
cation with published luminosity values in fields relatively 
"crowded" with Hll regions by Feinstein (1997) and Rozas 
et al. (2006). Velocity dispersion values and its errors are 
taken directly from Tables 2 and 4. The luminosity errors 
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have been estimated taking into account the errors in the 
amphtude (A) and the FWHM in the component fitting 
(F=1.0645*A*FWHM), and the distance errors. As a ref- 
erence value we have plotted a few Giant Hll regions from 
Bosch et al. (2002) together with their linear fit to their 
"young" GHll regions. 

Inspection Figure 11 reveals that the distribution of 
points in the \og[L) — log(cri„t) plane are strongly depen- 
dent on the components derived from the profile fitting. In- 
dividual components, labelled A and B (when applicable) 
as in their respective tables, and global values are those 
derived from a single Gaussian fitting to the line profile. 
As expected, individual components have smaller fiuxes and 
velocity dispersions than the global profile and points are 
therefore shifted in the diagram. The incidence of this analy- 
sis has, however, a different outcome for regions in NGC 6070 
and NGC 7479. In NGC 6070 the values derived from the 
global profile seem to lie, within observational errors, in the 
expected location if they follow the relation expected for viri- 
alised systems. When the individual components are plotted 
in the same diagram, only NGC 6070 II lies close to the re- 
gression, and NGC 6070 1 and IV seem to be too luminous 
for their measured velocity dispersion. On the other hand, 
two out of three regions in NGC 7479 lie well away this re- 
gression if parameters are derived from a single Gaussian, 
but four out of five individual components show a very good 
agreement when plotted in the same diagram. 

The interpretation of this behaviour is not straightfor- 
ward. Definitely, the presence of more than one Gaussian 
component rules the final position of the Hll regions in the 
log(I/) — log(cr) plane. Our high resolution spectroscopic data 



has allowed us to disentangle multiple kinematic compo- 
nents, which refiects in a big improvement in the x axis 
of the plot. But limitations in spatial resolution prevent us 
from identifying multiple components, if present, which do 
not become evident in the spectral signature. This has been 
a problem in previous attempts to derive relations in the 
log(I/) — log(cr) as discussed by Bosch et al. (2002) when 
they were able to obtain photometry for individual knots 
within Giant Hll regions in MlOl. 



4 SUMMARY AND CONCLUSIONS 

From new high resolution spectra of the Hll re- 
gions NGC7479I, NGC7479II, NGC7479III, NGC6070I, 
NGC6070 II and NGC6070 IV obtained at the 100-inch du 
Pont Telescope, LCO, we have confirmed the giant nature 
of all these regions. We have also found that all of them 
show evidence of wing broadening evident mainly in the Ha 
line and confirmed in other emission lines. In five of them 
we are able to fit a broad component which explains the 
integral profile wings. Only in NGC7479 II can we fit the 
integrated profile wings with two narrow components sym- 
metrically shifted in velocity with respect to the average of 
the two intense components. 

NGC7479 1 and NGC7479 II reveal the presence of two 
separated kinematic components with relatively narrow pro- 
files in all analysed emission lines. For NGC7479 1 region we 
can also generate 2D velocity images in Ha, [Nn]A6584A 
and [Sn]A6717A lines. These images show that the two nar- 
row components are spatially resolved and the kinematical 
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Table 4. Results of Gaussian profiles fitting to the observed emission lines in NGC 6070. Each emission line is identified by its ion 
laboratory wavelength and ion name in columns 1 and 2. According to the different fits performed on each line, column 3 identifies each 
"narrow component" (A and B, where applicable), a broad component and/or blue and red wings (b wing and r wing respectively). Radial 
velocities (Vr) and intrinsic velocity dispersions {(Tint) together with their respective errors arc expressed in kms~-^. The intrinsic velocity 
dispersions are corrected for the instrumental and thermal widths. Emission measures (EM) are included and shown as a percentage of 
the component flux relative to the total EM of the region. Upper panel: du Pont data and lover panel: Mike data. 
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information helps to split the knot in the velocity-distance 
plane. There is also evidence of the presence of a broad com- 
ponent in all the narrower extracted sections, at least in the 
Ha emission line. For NGC7479 II region we can also fit two 
relatively narrow Gaussian components. In both regions the 
radial velocity average value of our components A and B de- 
rived from the Ha emission line matches the global velocity 
in the galaxy when we project our value over the isoveloc- 
ity contours in Ha velocity field map. A similar comparison, 
using the HI map yields a value much closer to the value 
measured for component B. Then, we could only suggest 
that component A shows an odd kinematic behaviour. New 
observations with better spatial and spectral resolution are 
needed to clarify this point. 



In the rest of the studied Hll regions, we find one nar- 
row Gaussian component together with an underlying broad 
component, with no evidence of multiple narrow contribu- 
tions. The radial velocity average derived for the narrower 
component corresponds to the one expected from galaxy ro- 
tation curves, within the observational errors, in all of these 
cases. 

In most cases, the velocity dispersions of the broad com- 
ponents for these HII regions derived by optimal Gaussian 
fit are in agreement with Hagele et al. (2009, 2010) for cir- 
cumnuclear star-forming regions with atnt in H/3 line ~34 
to 65 kms~^, and those derived by Chu & Kennicutt (1994) 
and Melnick et al. (1999) for 30 Doradus nebula with aint 
in Ha line ~45 kms~^. 
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Figure 8. Ngauss fits witii two Gaussian components in tlie NGC6070I emission line profiles, du Pont data on the left side and MIKE 
data on the right side. In order from top to bottom panels: [Olll]5007A, H«, [Nll]6584Aand [Sll]6717A. 
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The estimated offsets between the narrow and broad 
components in NGC 7479 III, NGC 6070 II, and NGC 6070 
IV are in complete agreement with the values found in 
NGC 2903 and NGC 3310 by Hagele et al. (2009, 2010) re- 
spectively, which are between -25 and 35 kms^^. 



The distribution of the regions in the log(L) — log(cr) 
plane are strongly dependent on the components derived 
from the profile fitting. Individual components have smaller 
fluxes and velocity dispersions than the global profile and 
points are therefore shifted in the diagram. In NGC 6070 
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Figure 11. log(L) — log(cr) relation for our HII regions. Luminosi- 
ties and velocity dispersions are derived from our spectrophoto- 
metric data. The plot includes results from individual components 
labelled A and B (where applicable) together with velocity dis- 
persions for the overall profile (labelled as "global") which were 
derived from a single Gaussian fitting to the line profile (color 
solid error bars). The narrow components in NGC6070 were la- 
belled A, for which we only plot the MIKE measured data. A 
few Giant Hll regions from Bosch et al. (2002) (blue dashed er- 
ror bars) together with their linear fit to their "young" Giant Hll 
regions are plotted as a reference value. The luminosities are not 
corrected for extinction. 
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the values derived from the global profile seem to lie, within 
observational errors, in the expected location if they follow 
the relation expected for virialised systems. When the in- 
dividual components are plotted in the same diagram, only 
NGC 6070 II lies close to the regression, and NGC 6070 1 and 
IV seem to be too luminous for their measured velocity dis- 
persion. On the other hand, two out of three regions in 
NGC 7479 lie well away this regression if parameters are 
derived from a single Gaussian, but four out of five individ- 
ual components show a very good agreement when plotted 
in the same diagram. Definitely, the presence of more than 
one Gaussian component rules the final position of the Hll 
regions in the log(L) — log(a) plane. 
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Figure 9. NgauSS fits with two Gaussian components in the NGC6070II emission line profiles, du Pont data on the left side and MIKE 
data on the right side. In order from top to bottom panels: H/3, Ha, [Nll]6584A, and [Sll]673lA. 
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Figure 10. NgauSS fits with two Gaussian components in tiie NGC6070IV emission line profiles, du Pont data on the left side and 
MIKE data on the right side. In order from top to bottom panels: [Olll]5007A, Ha, [Nll]6584Aand [Sll]6717A. 
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